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Summary 


This report consists of two major parts: 1) the approach to develop morphing wing weight equations, and 2) the 
approach to size morphing aircraft. Combined, these techniques allow the morphing aircraft to be sized with 
estimates of the moiphing wing weight that are more credible than estimates currently available; aircraft sizing 
results prior to this study incorporated morphing wing weight estimates based on general heuristics for fixed-wing 
flaps (a comparable “morphing” component) but, in general, these results were unsubstantiated. This report will 
show that the method of morphing wing weight prediction does, in fact, drive the aircraft sizing code to different 
results and that accurate morphing wing weight estimates are essential to credible aircraft sizing results. 

The first portion of this report describes a method to develop morphing wing weight equations. The basic idea 
behind this approach is that a morphing wing weight database is necessary to develop a parametric wing weight 
predictor similar to those currently used for fixed-geometry aircraft components. Design of Experiments (DOE) 
provides a means with which to define a database of morphing wings with various planform dimensions, extents of 
shape variation due to morphing, and loading conditions. For each wing prescribed in this database, representative 
finite element models will then be created and structurally sized using an optimization routine developed for 
morphing wing components. The optimized weight of each of these FEMs, along with the parameters describing the 
wing, will be used to determine coefficients and exponents for a basis equation. The resulting equation becomes the 
morphing wing weight predictor. The means by which to implement this process using idealized beam model 
representations and intermediate complexity finite element models will be described. 

The second portion of this report describes the methods used to size a morphing aircraft; the major aircraft 
parameters and “optimal” wing shapes will be determined such that the aircraft gross weight is minimized for a 
given design mission. In this sense, the aircraft sizing process determines the basic wing dimensions and morphing 
capabilities necessary to most efficiently enable the aircraft design mission. The implemented technique will use the 
idea of a morphing wing “template” that restricts the wing shape variation during the course of the design mission to 
physically realizable mechanization strategies (an admissible substructure and actuation scheme capable of 
producing the changes in wing shape). The integration of the developed morphing wing weight equation into the 
template -based aircraft sizing routine will be shown to be a trivial matter in this case. 
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Nomenclature 


AR 

N 

S 

T/W 

W/S 

m 

n 

n z 

t/c 

{a} 

W 

{X} 

{X a .J 

Ab 

Ac 

AA 

X 

A 


Wing aspect ratio 
Number of DOE trials 
Wing area 

Aircraft thrust to weight ratio 
Wing loading 

Number of mission legs in the aircraft design mission 

Number of DOE system parameters 

Maximum vertical load factor 

Wing thickness to chord ratio 

Set of fitting parameters in the basis equation 

Structural design variables 

Set of morphing wing parameters 

Set of aircraft sizing design variables 

Change in wing span 

Change in wing chord 

Change in quarter-chord wing sweep 

Wing taper 

Quarter-chord wing sweep 


I. Introduction 

M ORPHING aircraft, as defined in this report, are aircraft utilizing wings that have the capability to drastically 
change planform shape during flight - perhaps a 200% change in aspect ratio, 50% change in wing area, and 
a 20 degree change in wing sweep. 1 This type of design might be incorporated to enhance various operational 
capabilities of the aircraft, reduce the aircraft’s required takeoff gross weight, and/or enable an aircraft to fly a 
design mission that a fixed-wing aircraft could not.. Most recently, two flight-traceable morphing wing concepts 
were developed through Phase II of DARPA’s Morphing Aircraft Structures (MAS) program. 2 The first, a “folding 
wing” concept, was developed by Lockheed Martin and enables variations of span length, aspect ratio, and effective 
sweep angle. The second, a variable sweep / variable root chord concept, was developed by NextGen Aeronautics 
and enables direct variations in root chord length and sweep angle; indirectly varying the planform area and aspect 
ratio. Both of these are illustrated in the following figure. 

Figure 1 Recently Developed Morphing Wing Concepts (Left - Lockheed Martin’s Folding Wing Concept, 
Right - NextGen Aeronautics Variable Sweep / Variable Root Chord Concept) 

While Phase II of the MAS program concentrated on development and testing of scaled wind tunnel models to 
determine system feasibility in flight equivalent environments (wind tunnel experiments at the TDT at NASA, 
Langly), Phase III will pursue the development of these concepts as flight demonstrational vehicles. 3 

Clearly, application of morphing wing technologies will require the concurrent development of design and 
optimization strategies at the aircraft level to expedite overall development of these systems. Of major importance 
is the development of a robust morphing aircraft sizing code to be used during conceptual design tasks. Such a tool 
would enable studies of the operational benefits at the aircraft system level and provide a methodological basis for 
future morphing aircraft sizing codes. The major drawback of such an effort, however, is that accurate weight 
predictions are required for all major components, including the morphing wing. Little information is currently 
available to substantiate such a wing weight prediction and thus, the sizing results would be unfounded. However, 
this report illustrates various means to implement a robust morphing aircraft sizing code (Part III) and develop 
morphing wing weight predictions that are more credible than currently available (Part II). The entire process is 
then demonstrated using a variable sweep / variable chord concept similar to the NextGen concept. 


3 


II. Developing Finite Element Based Wing Weight Predictors 

The most common approach to developing fixed-wing weight predictors is centered on the idea that a large 
database of wing weights, their associated geometry, and maximum loading conditions already exists through 
fielded aircraft. By selecting a subset of these data with common features (e.g. wings having aspect ratio and taper 
within some range), various analytical methods can be employed to predict flight-ready wing weight as a function of 
the wing parameters (e.g. AR, W/S, n z , t/c, etc.). Reference 4 provides an early example of this technique (circa 
1968) in which fighter aircraft wing weight predictors were developed from a small set of aircraft (~15) using a least 
squares regression approach and basis functions similar to those developed by F.R. Shanley 5 in 1960. 

Direct application of the above approach is not possible for morphing wing components because an adequate set 
of fielded aircraft data is nonexistent. However, the approach used here is to develop a surrogate wing weight 
database, and proceed with the equation development in the usual way. This sort of approach was successfully used 
in the design of a high speed civil transport aircraft (fixed-wing). 6 Using Design of Experiments (DOE) methods to 
define a set of morphing wings with various shapes, representative finite element models are developed for each 
wing in the database and then sized to give a corresponding weight estimate. These data are then “best fit” to an 
appropriate basis equation using a least squares regression technique resulting in the morphing wing weight 
equation. The various aspects of this process are described in more detail in the following sections. 

A. Design of Experiments 

In this research, a morphing wing will be defined by a number of parameters, {X}, that include the basis wing 
geometry (also referred to as the reference wing geometry), extent of morphing capabilities, and parameters 
describing the structural loading conditions. More specifically, {X} is defined by the following relation: 

{X} e { AR, S, t/c, A, A} u { Ab, Ac, AA}u{W/S, n z } (1) 

where the sets have been separated by category. Although {X} must contain the basis wing geometry and loading 
conditions, the set of morphing parameters contained in {X} depends on the particular morphing wing concept being 
implemented. For example, a conventional variable sweep wing would be represented by the morphing parameter 
set {AA} because the morphing mechanism directly varies sweep. Note that although a variation in sweep changes 
the chord distribution and span length measured with respect to the aircraft body axis, Ab and Ac are not included in 
the set because the morphing mechanism has not physically changed the span and chord lengths measured in the 
wing’s reference frame. 

The Design of Experiments technique enables the mapping of a continuous design space to a finite set of points 
within the design space. Various Designs of Experiments exist, and the particular type used governs the number of 
points in the mapped set and their distribution / position within the design space. Table 1 illustrates three possible 
DOE types and the number of “experiments” (points within the mapped design space) as a function of the number of 
independent parameters to be varied (the number of dimensions making up the design space). For example, if a 
moiphing wing was described by a set of 7 parameters, a full factorial design would prescribe 2187 wing shapes 
within the bounded design space at which to evaluate the wing weight, while a face-centered central composite 
design would prescribe 77. The table entry for the saturated D-optimal design shows two values for the number of 
prescribed experiments. This results from the fact that the number of design points specified is equal to the number 
of unknowns in the basis equation used during the regression. As will be shown in Section D, two such basis 
equations are considered and have different numbers of unknowns. For any of these table entries, the left number is 
representative of the “conventional” basis equation, while the right represents the full quadratic response surface. 



Number of Parameters 

DOE Type 

6 

7 

8 

n 

Full Factorial Design 

729 


6561 

3" 

Face-Centered CCD 

77 

143 

273 

2"+2n+l 

Saturated D-optimal Design 

7/28 


9/45 

(n+1) / (n+l)(n+2)/2 


Table 1 Number of Experiments for Various DOE Types and Number of System Variables 


Typically, the more data available from the design space will reduce the uncertainty of wing weight variations 
within the design space, but will clearly increase the time needed to develop the weight database. Therefore, the 
goal of selecting a DOE type is to maximize the quality of the regression (accurately determine the wing weight 
trends with respect to variation of parameters) while minimizing the number of points used to formulate the weight 
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equation. Unfortunately, the quality of regression typically can not be predicted a priori and the convergence of the 
weight equation must be determined through trial and error. However, in Reference 7, this author finds that a 
saturated D-optimal design set gives a good representation of trends within the design space, while a face-centered 
central composite design yields equations with coefficients of correlation on the order of 0.98. Note that this 
heuristic was based on the “conventional” basis equation (see Section D) using beam model representations of the 
morphing wings. 


B. Structural Optimization Strategies 

1 . Simultaneous Analysis Approach 

Although methods to carry out the structural sizing process are outlined in greater detail in previous work (see 
Reference 8), the simultaneous analysis approach has been used in this investigation and will be revisited here for 
completeness. The aircraft is assumed to operate with the morphing wing in distinct configurations (shapes) for 
major portions of the design mission. For example, the wing might operate in a high sweep configuration during a 
dash leg of the design mission. In this approach, the structural sizing routine follows a conventional strategy for a 
fixed-wing structure in which a set of design load conditions (e.g. symmetric and / or asymmetric maneuvering 
conditions, gust loads, impact loads, etc.) are considered simultaneously during the sizing process and an “optimal” 
structural weight is determined such that all imposed design constraints are satisfied. In the case of a morphing 
wing, the critical load cases which tend to induce maximum constraint conditions must be considered 
simultaneously for all wing configurations to determine the optimal (or nearly optimal) load paths. Figure 2 
illustrates this methodology with greater 
detail. Here, {x} represents the structural 
design variables (i.e. spar thicknesses, 
spar cap areas, number of layers in a 
composite laminate, and/or skin 
thicknesses) and {g}j represents the 
design constraints (i.e. maximum stress, 
buckling, maximum deflections, etc.) for 
each wing configuration i. This type of 
sizing allows for meaningful aeroelastic 
constraints to be imposed, such as flutter 
tendencies, which might be difficult to 
enforce in other structural sizing 
approaches; see, for example, the 
“aggregate” sizing approach in 
Reference 8. 

2. Simultaneous vs. Sequential Optimization Strategies 

Although the simultaneous analysis structural sizing approach is robust and follows a traditional structural sizing 

approach, its implementation is not readily implemented with common FEA toolsets. A secondary strategy, as 

presented by Lockheed Martin, is to size the wing using a sequential approach. 9 In this case, the wing is sized 
subject to a set of load cases in a load dominant configuration using a common FEA/optimization environment. The 
sized structural dimensions are then used as lower design variable bounds in subsequent structural optimizations 

with the wing in non-dominant load configurations. In the case of the Lockheed morphing wing design, the idea 
seems intuitively clear because the morphing is relatively simple (having only two configurations) and the 
substructure is fixed; the problem is simple enough to rationalize the use of the sequential optimization. However, 
in a more complicated structure - one which has a spatially variant substructure - with multiple configurations and 
load conditions, it is possible for a non-dominant load configuration to drive an optimal load path due to the 
arrangement of the substructure. In fact, for a simple morphing wing model, this author finds that the sequential-like 

optimization approach tends to over-predict the structural weight by an average of 13.5% when this sort of situation 
occurs. 2 * * * * 7 



Figure 2. The Simultaneous Analysis Sizing Approach 
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C. Finite Element Modeling 

The level of detail in the finite element model (FEM) representation of the morphing wing will play a pivotal roll 
in developing accurate weight predictors. Two FEMs are considered here: 1) a simplified beam model representing 
a low fidelity model and 2) an intermediate complexity model. The beam model representation allows rapid 
development and validation of the techniques (optimization and regression) needed to create the weight equation, 
can quickly perform trade studies for a given morphing wing concept, and provides a basis for comparison with 
more complex models. The intermediate complexity FEM will more accurately model the wing’s substructure, 
provide a more accurate prediction of the fully designed wing weight (leading to a more reliable weight predictor), 
and indicate the system aspects that are most sensitive to modeling fidelity when compared to the beam model 
representations. 


1. Structural Modeling 

a) Beam Model Representation 

The beam model representation of the morphing wing provides a means to represent and size for the required 
bending material of the wing; this accounts for a majority of the weight when considering a fixed-wing structure 5 . 
To approximate the morphing wing structure, a beam is defined along the quarter-chord line of the wing in its 
reference configuration (previously referred to as the “basis” geometry). This beam is discretized into a number of 
segments such that the taper of the wing can be accounted for; the thickness of any segment is the average of the 
wing thickness at the element’s endpoints. Each discretized segment of the beam is treated as a beam element in the 
FEM. The cross-section of each element is rectangular with the height specified by the wing thickness and the 
width being a design variable during the structural optimization routine. 

The quarter-chord beam is defined in every model. Flowever, this might not be sufficient to accurately describe 
the envisioned morphing mechanism. For example, a secondary beam might be modeled to provide a mechanical 
linkage to move and support the main beam. These secondary beam components are also modeled, but their 
geometry and discretization are dependent on the specific wing concept being implemented; References 7, 8, and 10 
provide several examples of this. 

Aerodynamic loads are predicted for a symmetric pull-up maneuver with specified design load factor (the design 
load factor is defined as the limit or operational load factor multiplied by the safety factor of 1.5) via the vortex 
lattice method (VLM) presented in Ref. 11. For any given element making up the main beam, the resultant vertical 
load predicted over the portion of the wing span bounded by the element’s span domain is reduced to equivalent 
nodal loads at the element’s endpoints; these loads are applied to the element neglecting the angle of attack of the 
wing (no horizontal loads are applied to the element). Furthermore, because this simple structural model does not 
accurately account for a real wing’s torsional stiffness, chordwise pressure variations and the resulting structural 
torsion loads are neglected. 


b) Intermediate Complexity Model - Morphing Wing Templates 

The intermediate complexity models introduce a more detailed representation of the wing’s substructure by 
modeling a box-beam or shell type structure capable of spatial rearrangement to facilitate moiphing. In this case, 
the model and subsequent analyses are not as restricted as in the beam model case. For example, the critical load set 
is not restricted to the single symmetric load case, but can include a complete set of symmetric and antisymmetric 
maneuvers. Furthermore, because the model now more accurately reflects the true structure of a morphing wing, the 
sized structure should give a better representation of the morphing wing’s true weight relative to the representative 
beam model; provided that the main weight drivers of the wing have been accounted for in the FEM. 

Development of an intermediate complexity model (ICM) currently revolves around the idea of a morphing wing 
“template”. In the most basic sense, a wing template defines a structural mechanization strategy that enables 
specific shape changing capabilities of the wing. For example, the variable sweep wing found on the F-14 could be 
referred to as a rigid-body sweep template. This template would describe a fixed inboard wing structure, a pivot 
joint, and the variable sweep outboard structure. The driving idea behind this “template” -based approach will 
become apparent in Section III when morphing aircraft sizing strategies are discussed. 
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Development of an ICM is currently facilitated using a series of Matlab-based scripts. The overall process flow 
is illustrated in the following figure: 

Matlab Subroutine 


Initialize 


Template 


Discretize 


Extrude 


Morph 


Write FEM 


Figure 3 Process Flow for the Development of a Template-based FEM 

The overall process requires only the definition of the morphing wing template and the discretization scheme while 
the other routines are fully autonomous. Also note that this scheme enables the rapid development of a morphing 
wing FEM as a function of some baseline wing geometry and the extent of actuation. This is important because the 
DOE task requires the rapid evaluation of many wings of varying sizes and morphing capabilities; a nontrivial task 
if one were to manually develop each model. Currently, several intermediate-complexity FEMs have been 
developed with various templates, but have not yet been implemented in a structural sizing process. Appendix A 
illustrates several examples of these. 



Description 


Define basic planform geometry based on 
baseline wing geometry and template type 


Operate on basic geometry: discretize 
planform shape to static and morphing 
sections, define spar / rib layout 


Segment spar / rib layout mesh as desired; 
these define element bounds in 3D FEM 


Extrude the segments to the upper and 
lower wing surfaces to define nodal 
coordinates and elements 


Operate on morphing wing sections 
{stretch, shear, rotate, translate} 


Output FEA input decks for all wing 
configurations 


2. Major Structural Joints 


An important aspect of the morphing wing is that 
rearrangement of the substructure will require structural joints 
to facilitate the overall shape change. For example, the 
variable sweep capabilities of the F-14 require the use of a 
substantial pivot pin to connect the wing and fuselage 
structures. In this work, all pivot joints are idealized to the 
model illustrated in Figure 4. The design task is then limited to 
preventing failure due to shear in the pivot pin and the sized 
pin will represent the entirety of the pivot joint weight. The pin 
is assumed to be a solid cylinder (grey in the figure) made of 

Aluminum 2024-T6. The bending loads from the outboard wing section (green) are reduced to a force couple with 
magnitude P and moment arm 3h/4. Flere, h is specified by the thickness of the wing at the pivot joint. The pin 
diameter is then sized to prevent plastic deformation due to the maximum shear stresses on the circular cross- 
section. Reference 7 provides a more detailed representation of this model. 



FBD of Pin 


^r F 

0.75h 


Figure 4 Structural Pin Model 


3. Actuator Weight Prediction 

A further important aspect of a morphing wing is the actuating mechanism that physically rearranges the wing 
structure. Currently, this mechanism is represented by an idealized hydraulic actuator modeled as a simple pressure 
vessel (the cylinder) and a solid circular column (the ram). The actuator loads are determined by adding the effects 
of the loads required to rearrange the aerodynamically loaded structure from one configuration to another (which is 
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related to changes in structural strain energy), and the loads required to overcome friction induced by moving 
structural components (e.g. rotation of the loaded pivot joint). These loads are then used to size the ram subject to 
strength and buckling constraints while the cylinder is sized subject to strength constraints using a nominal fluid 
pressure value. Further details regarding the modeling of these mechanisms can be found in the Appendix of 
Reference 7, with various case studies presented in Reference 8. 

D. Weight Equation Basis Functions 

The basis function is a generic equation form with unknown parameters that must be “best fit'’ to some given 
data set. In this work, two such basis functions will be considered when formulating the wing weight equation: The 
full quadratic polynomial and the conventional equation form (a product of terms raised to unknown constant 
powers). Table 2 illustrates both of these equations. 


Equation Type 

Functional Form 

Number of Unknowns 

Full Quadratic 
Equation Form 

>r„, ({*,„}) = +Z2>AA 

i=\ i=l j=i 

(n + l)(n + 2) 
2 

Conventional 
Wing Weight 
Equation Form 

i= 1 

(» + 1 ) 


Table 2. Functional Forms of the Wins Weisht Basis Eauations 


Here, X-, represents one of the n system parameters within the DOE set, {X}, and a, are the unknown constants that 
must be fit to the data set. The conventional equation form arises from an analytical development of the wing’s 
bending material weight as a function of the DOE parameters, {X} (See Reference 7, Appendix B). This author has 
found that the full quadratic equation form does not accurately capture / reflect the trends within the design space 
and leads to a poor correlation when fit to the known data (regardless of the number of trials in the DOE). On the 
contrary, the conventional equation form tends to adequately model the weight trends, even for a saturated DOE trial 
set, and is the suggested basis equation for this type of application. 

E. Regression Technique 

A least squares regression is used to fit the data to the basis function. More specifically, given values for the 
fitting parameters, {a}, the approximated wing weight for each wing shape, [X j-,, from the DOE database is 
evaluated using the wing weight predictor. The difference between the actual wing weight from the FEA 
optimization and the approximated weight from the wing weight equation is calculated and squared; this is referred 
to as the squared error. The squared error is calculated for each wing in the DOE database and then summed to give 
the total error, or equivalently, the sum of squares. The optimal values for {a} are those which minimize the sum of 
squares and thus “best fits” the data to the basis equation. 

F. Process Implementation 

1. Beam Model Representation 

The beam model representation was the first considered by this author. In order to expedite development of the 
aforementioned techniques and quickly perform various trade studies, the entire process was implemented in a 
common environment. In this case, all tools were developed, integrated, and executed within Matlab 12 , including 
the FEA code. The entire process layout is illustrated in Appendix B of this document. 


2. Intermediate Complexity Model 

The intermediate complexity model requires the use of a more advanced FEA environment like ASTROS 13 or 
NASTRAN 14 . Therefore, the overall process is more complex - requiring the integration of potentially multiple 
software environments (DOE, FEA, optimization, and regression toolsets) - yet more realistic to implementation 
outside the academic environment. Currently, iSIGHT 15 integrates the various process components necessary to 




implement the “simultaneous analysis” structural optimization routine. This program implicitly provides a means to 
integrate various toolsets by: Controlling the main design variables of the problem (e.g. evaluating a DOE trial set 
or carrying out an optimization task), outsourcing tasks to other software packages, and autonomously handling the 
file I/O between itself and the various programs. In essence, this program automates trade studies and portions of 
the design cycle. Appendix C illustrates the overall process flow for this model. 

G. Case Study - Variable Sweep / Variable Chord Concept 

The following section will briefly describe a case study in which a wing weight equation is developed for a 
concept similar to NextGen Aeronautics morphing wing concept. References 7 and 8 should be consulted for 
further modeling details as this study will parallel those found in the references. Any differences will be clearly 
defined and the resulting equation will be utilized in the aircraft sizing case study of Part III of this report. 


1. Definition of the Wing Geometry’ Parameters, {X} 

The variable span / variable chord wing concept is characterized by the following parameter set: 

{X} e {AR, S, t/c, A, A}u{Ac root , AA}u{W/S} (2) 

where the variable n z is set to a limit load factor of 1.5 for all wing configurations to reduce the number of DOE 
parameters. 

2. DOE Parameters and Their Bounds 

Because seven parameters are being used to describe this wing concept, 273 DOE experiments will be required 
to implement the FC-CCD set (see Table 1). The variable bounds in the DOE set are defined as: 



Bounds 


Parameter 

Lower 

Nominal 

Upper 

Unit 

AR 

5 

10 

15 

[1 

S 

200 

300 

400 

[ft A 21 

t/c 

0.06 

0.09 

0.12 

[1 

X 

0.3 

0.4 

0.5 

[1 

A 

5 

15 

25 

[deal 

Ac root 

0.5 

0.7 

0.9 

[1 

AA 

5 

17.5 

30 

[deg] 

W/S 

50 

100 

150 

[lb/ft A 21 


Table 3 DOE Parameters and Their Bounds 


3. De finition of the Wing Con figurations and Structural Layouts 

Figure 5 illustrates the half-span planform geometry (red) and structural elements (blue) in this model using the 
nominal variable bounds defined in Table 3. Four distinct configurations are defined and are intuitively associated 
with particular flight regimes as indicated by the plot titles. Also note that the forward spar (beam along the quarter- 
chord line) is discretized using 5 elements while the aft spar (beam controlling the position of the forward spar) is 
modeled using a single element. Each element thickness is a design variable during the structural optimization 
routine. 
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0 


Ugh Lift Config 


Dash/Cruise Oonfig 



S : 127.21 [ft 2 ] 

AR: 4.2 

s ' / ' fiep c/ 4 : 35.° [deg] 
Taper : 0.40 



4. DOE Results 

The mean and standard deviations of wing weights computed for the DOE trials are summarized in the following 
table. Across the broad range of factors, the standard deviation is large compared to the mean value of weight. 


Sample Metric 

Value 

Unit 

Mean Weight 

2870 

[lb] 

Std. Deviation 

3310 

[lb] 


Table 4 DOE Results 


5. Regression Results 

The following Figure illustrates the regression results for the conventional equation form and the full quadratic 
equation. Here, the structurally sized weight resulting from the FEA optimization is plotted against the predicted 
weight of the wing using the regressed wing weight equation. The heuristic proved applicable in this case with the 
conventional equation form fitting the data quite well (a coefficient of correlation of 0.9959). The quadratic 
equation was unable to predict the data well near the extreme weights (a coefficient of correlation of 0.9783). The 
coefficient of correlation is a measure of the linearity of the data in the following graphs. 
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RJI Quadratic Equation 


Actual VNfeight [lb] 



Figure 6 Comparison of morphing wing weight equations; conventional wing weight equation form (left) and 

quadratic equation form (right). 


The resulting conventional-basis equation is: 

tO.4516 t * — \ 1.3520 


W wilE =(0.0352- W c ). 


S“ 5i6 (AR) .(1 + A) 


0.7850 


1 + 


AA 

90 


•(1 + AC max) 


0.0252 


(l + cosA ) 4 ' 8970 -(lOO-f/c) 


0.8870 


f w/s ^ 00045 
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The units for these variables are the same as in Table 3 and the resulting wing weight is in pounds. 


(3) 


TIT. Preliminary Morphing Aircraft Sizing 


A. Overview 

Conceptual sizing methods for fixed-geometry aircraft predict the designed weight of an aircraft to perform a 
given mission based on the values of major aircraft parameters. Six such parameters include wing loading ( W/S ), 
aspect ratio (AR), wing sweep (A), thickness-to-chord ratio (t/c). wing taper ratio (1), and thrust-to-weight ratio 
( T/W ); these particular variables are referred to as the “basic six” because they typically have the greatest impact on 
aircraft size, weight, and overall flight performance 16 . The conceptual design task, then, is to determine the values 
of these design variables such that the “sized” aircraft meets all performance requirements at a minimum takeoff 
gross weight and/or estimated cost; this is typically facilitated with trade studies or formal optimization techniques 
(gradient-based or global search algorithms). 

For a morphing-wing aircraft, the conceptual design task substantially increases in complexity; in addition to the 
basic six aircraft parameters, the “optimal” variation in wing geometry throughout the design mission must also be 
determined. The approach used to augment the fixed-wing aircraft sizing task to account for the morphing wing 
capabilities is the topic of this section. 

B. Problem Statement 

Although the morphing capabilities will increase the number of independent variables during the aircraft sizing 
routine, the overall problem statement remains the same: 
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Minimize: 


(Takeoff Gross Weight) 


(4) 


W G 


Variables* 111 = 1 w/s > T/w } u { AR , S, t/c, X, A} u {Morphing Parameters} 

Mission Performance Constraints 


Subject to: 


(5) 


. Takeoff Parameter < 1 15 

• Cruise > Mach 0.7 

. Landing Distance < 5000 [ft] 

• Excesss Power > 0 for each mission leg 

Variable Bounds 

• lbj < {X a .c.j} < ubj 


(6) 


Note that the performance constraints as defined above are neither necessary nor exhaustive, but provide an 
illustration of the problem statement. The set of wing morphing parameters are the discussion of the next section. 


C. Independent Variable Sets 


1 . ‘ 'Photomorphing ” Approach 


In past research efforts, the traditional sizing approach is modified to incorporate morphing wing capabilities in 
an approach referred to as “Morphing as an Independent Variable”. 171819 Essentially, a mechanization strategy is 
not defined for the moiphing capabilities and the wing geometry is allowed to freely “morph” into different shapes 
and sizes for each phase of a defined mission profile. This concept of wing morphing has been referred to as 
“photomorphing”, after the popular graphics software program named Photoshop, in which images can be altered 
without regard for the physical practicality of the resulting image. This approach requires that the “basic six” design 
parameters related to the wing geometry be expanded to individual variable sets for each leg of the design mission. 
Therefore, the independent variable set within the general aircraft sizing problem statement becomes: 


Independent 

Variables: 


( x a.,} = {W/S, T/W}u{AR ; , S i? t/c„ A„ A, 


i = 1 :m 


(7) 


where m is the number of mission legs. Variation of all the wing parameters for each mission leg is not a necessary 
condition for this problem setup. For example, the following definition of independent variables is also admissible: 


Independent 

Variables: 


{X, C .MW/S, T/W}u{A}u{AR i; S„ t/ Cj , T,} i 


l:m 


( 8 ) 


In this case, the quarter-chord wing sweep is still a variable in the overall problem, but it is constant for each mission 
leg; the other four wing parameters still vary for each leg of the aircraft design mission. 

As highlighted in References 18 and 19, a drawback to this methodology is that the independence of the 
geometric design variables does not follow any particular mechanization strategy for a morphing wing; solutions to 
this problem (indicated morphing capabilities) might not be physically implemented. A further concern is that the 
morphing wing weight predictions rely upon wing weight equations for fixed-geometry aircraft, with a general 
penalty added to account for structural actuation (See Section 1). While it is unclear how this prediction would 
compare with an actual structure, it is clear that the accuracy of the prediction will drive the quality of the sizing 
results. 


2. Template-based Approach 
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A specific morphing wing concept will most likely be restricted in its shape changing capabilities because of the 
internal mechanization of the wing’s substructure. The template -based sizing approach builds on the photomoiphing 
technique by incorporating this idea. The basic concept underlying this approach is the definition of a reference 
wing geometry. This geometry is described by a set of parameters that uniquely define the overall geometry of the 
wing. Continuing with the standard set of variables defining the wing, the reference wing geometry will be defined 
as follows: 

{X ref } = {AR, S, t/c, A, A} (9) 

The overall geometry of the morphing wing can now be defined by augmenting the reference geometry with 
parameters describing the shape variation attributed to morphing, {X morph } : 

{ X morp h }e{Ab, AC,AA} (10) 

where it is noted that the set can be defined using a single parameter, or any combination thereof because they are 
linearly independent. Here, Ab and Ac represent the percentage increase in span and chord length with respect to the 
corresponding reference geometry values. The AA term is a normalized parameter describing the change in quarter- 
chord sweep angle. Mathematically, these parameters could define a morphed wing geometry described by the 
following functions: 

b, = (l+Abj ) - b ref 
Ci = (1 +Ac, ) - c ref 

X i = A ref +(AAj) -(;r/2) on 


S; = /(bpCjjAj) 

AR ; =/(b i ,c i ,A i ) 

Other wing parameters, {t/c, k}, would be updated as a function of the parameter values similar to S; in the previous 
equation. Thus, the definition of a reference geometry parameter set, moiphing variable set, and update formulas to 
describe a “morphed” wing shape is synonymous with the definition of a “template” in Part II. C of this report. 
However, in this case, knowledge of the internal mechanization strategy of the substructure is unnecessary, but still 
reflected in the geometry update formulas (equation set 9). 

The template -based aircraft sizing approach requires that the morphing parameters be expanded to individual 
variable sets for each leg of the design mission. Therefore, the independent variable set within the general aircraft 
sizing problem statement becomes: 


Variables^ 111 l X -i s l W/S > T/W } u { x «}^{ x »^}, (12) 

{Xmorph},. e{Ab„ Ac,., AA,} (13) 

where i represents a leg of the aircraft design mission, and m defines the number of mission legs. 

As is currently defined, the following formulation for wing geometry in terms of a reference geometry and a set 
of morphing parameters could lead to a non-unique design space; that is, multiple sets of wing parameter values 
could lead to the same morphed wing shape, leading to an ill-conditioned aircraft sizing problem. Therefore, the 
following constraints are imposed to avoid this condition: 

m 

Abj > 0 f^Ab,=0 


m 

QAc,=0 (14) 

i=l 


Ac- > 0 


n 


AA, >0 


TJaAi = o 

i=l 
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These constraints ensure that the reference wing geometry represents the minimum morphing parameter value (i.e. 
zero if the equality constraints are satisfied) for at least one of the mission legs. Equivalently stated, the reference 
wing reflects the minimum span length, chord length, and sweep angle the wing will take on over the course of the 
mission. Notice that this formulation does not force us to impose the reference wing geometry at a given mission 
leg, but does allow the reference wing geometry to be an admissible shape at some point during the mission. 

D. Aircraft Weight Prediction 

The gross weight of the aircraft is predicted during the aircraft sizing routine via the following equation: 

W G = W Empty + W Fuel + W crew + W Fixed ( 1 5 ) 

where the aircraft empty weight is predicted using parametric equations for the major structural components (e.g. 
fuselage, wing, horizontal/vertical stabilizers, pylons, etc.), the fuel weight is determined by simulating the flight of 
the aircraft design mission via governing equations of flight and general heuristics for fuel consumption, the crew 
weight is determined by mission specifications, and the fixed weight accounts for any miscellaneous items that 
should be accounted for in the weight prediction. Notice that the only nontrivial weight estimate is that for the 
moiphing wing. Methods to approximate the morphing wing weight are the subject of the following sections. 


1. Wing Weight Prediction via “Photomorphing” 

In the case that the wing is allowed to change shape via the photomorphing approach, the wing weight is 
calculated using a fixed-wing weight equation and then adding a weight estimate for the extent of morphing 
capabilities. More specifically, the fixed-wing weight is calculated for each wing shape during the design mission, 
and the representative fixed-wing weight is the maximum of these. The weight estimate for the morphing 
capabilities is determined by calculating the maximum variation in planform area from one mission leg to the next 
and applying a weight factor 5 pounds per square foot to this change in wing area; this is described in further detail 
in Reference 17. 


2. Wing Weight Prediction via Template-based Morphing Wing Weight Predictor 

Given that the aircraft sizing process is setup via the template-based approach, a response surface based wing 
weight equation can be developed through the methods defined in Part II of this report. Direct application of this 
equation in the empty weight prediction is then a trivial task, requiring only the reference wing geometry and 
maximum morphing parameter values; both of which are included in the aircraft sizing design variable set. 


E. The Optimization Problem 

The aircraft sizing algorithm is implemented within Matlab using fmincon - a gradient-based optimization 
routine using an SQP algorithm. The problem statement is the same as displayed in Section B except for an 
additional constraint that forces the gross weight of the aircraft (a design variable) to be within one pound of the 
calculated aircraft gross weight (sum of empty weight, fuel weight, crew weight, and fixed weights). This is 
illustrates in further detail in Reference 17. 


F. Case Study 

The following case study illustrates results of the described methods using a variable sweep / variable chord 
moiphing wing concept similar to the NextGen Aeronautics concept. 


1. Mission Definition 

The following figure and table defines the aircraft design mission for this case study. Note that the wing shapes 
for the decent legs have no bearing on the sizing results; a fuel fraction heuristic is used to determine fuel weight and 
no range credit is allotted for this leg such that the wing shape has no effect on the aircraft gross weight. Therefore, 
design variables are not needed for these mission legs. 
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2. Definition of the Wing Template 



The variable sweep / variable chord concept is capable of direct variation of the quarter-chord sweep angle and 
root chord length. Indirectly, changes in planform area, aspect ratio, taper, and thickness-to-chord ratio are also 
realized. The following sets define the wing geometry and its morphing capabilities: 


Leg 

Description 

Details 

Constraints 

1 

Startup, Taxi, Takeoff 

Altitude 

101 ft 

TOP <= 115 

2 

Climb 

Altitude 

[0 30]x1000 ft 


3 

Cruise 

Altitude 

[30 30]x1000 ft Leg Distance: 1500 nm 

Cruise Speed >= Mach 0.7 

4 

Loiter 

Altitude 

[30 30]x1 000 ft Endurance: 4 hrs 


5 

Decent 

Altitude 

[30 5]x1000 ft 


6 

Dash 

Altitude 

[5 5]x1000ft Leg Distance: 250 nm 

Dash Speed >= Mach 0.8 

7 

Climb 

Altitude 

[5 301x1000 ft 


8 

Cruise 

Altitude 

[30 30]x1 000 ft 


9 

Decent 

Altitude 

[30 0]x1000 ft 


10 

Land, Taxi, Shutdown 

Altitude 

[0]ft 

Landing Distance <= 5500 [ft] 


Table 5 Mission Leg Details 


{X„ r } = {AR,S,t/c, 2, A) 

(16) 

{X„„ >1 }={AC„,AA} 

(17) 


Two variables are needed for each mission leg to describe the morphing parameters. Therefore, a total of 16 design 
variables are needed to define the shape variation for the entire mission (the wing shapes for the decent legs are not 
designed). Other modeling assumptions include: 

• The length of the quarter-chord line in the reference wing geometry is invariant 

• The thickness to chord ratio remains constant during morphing 

• The change in sweep angle (AA) refers to the change in quarter-chord sweep angle 

• The chord at the wing tip remains constant in length and is oriented parallel to the flow (this was 
also assumed during the development of the wing weight equation) 


3. Sizing Results 

a) Photomorphing Results 

In the case that the photomorphing wing weight estimate is used in conjunction with the variable sweep / 
variable span wing template in the aircraft sizing routine, the sized aircraft has the following characteristics: 
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Weights 



W G 

32489.5 

[lb] 

W Fuel 

14959.3 

[lb] 

W Empty 

17030.1 

[lb] 

W Wing 

1769.5 

[lb] 


Table 6 Weight Summary of Sized Aircraft using Photomorphing Wing Weight Estimate 

The wing geometries are illustrated in Figure 8. 


b) Variable Sweep / Variable Span Wing Weight Equation 


In this case, a morphing wing weight equation was developed for the variable sweep / variable span template 
using the methods from part II. The developed wing weight equation, which is a function of the reference wing 
geometry and the morphing parameters ({AA, Ac r00t }), has the following form: 


S“ 16 (AR)'“° (1 + A) 


W„, s = (0.0352 -W 0 ) 


\ 0.7850 


' AA ^°- 4524 


1 + 


90 


(l + AC max) 


0.0252 


(1 + cos A) 4 ' 89 ™ -(100-f/c) 


0.8870 


w/s 

10 


, 0.0045 


(18) 


where the following basic units are assumed (ft, lb, deg). Using this wing weight equation, in conjunction with the 
variable sweep / variable span wing template during the aircraft sizing routine, the resulting sized aircraft has the 
following characteristics: 


Weights 



W G 

26495.0 

[lb] 

W Fuel 

11662.4 

[lb] 

W Empty 

14332.6 

[lb] 

W_Wing 

1767.7 

[lb] 


Table 7 Weight Summary of Sized Aircraft using the Variable Sweep / Variable Span Wing 

Weight Equation 


The “optimal” wing geometries appear in Figure 8. 


c) Summary of Results 

The following figure illustrates the wing planform shapes for the resulting aircraft sized using the 
photomorphing wing weight prediction and the response surface based wing weight equation. 
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Figure 8 Superposed wing planform sizing results using the Variable Sweep / Variable Chord Template (Red - 
Photomorphing wing wt. prediction; Cyan - Response surface based wing weight equation) 


The wing weight prediction alone has a significant impact on the optimal shape of the wing throughout the 
design mission. This is an important result, because it suggests that an accurate wing weight prediction is necessary 
to determine the overall size of the wing and the required extent of morphing actuation during the preliminary 
design phase. Substantial deviation of the fully designed wing weight from the estimated wing weight during 
preliminary design could lead to an aircraft unable to meet the overall mission requirements. 


IV. Conclusions and Suggested Future Efforts 

The overall process of developing a morphing wing weight equation and directly incorporating the weight 
estimate into a template-based sizing code tends to be a viable methodology. The particular morphing wing weight 
estimation used during aircraft sizing was shown to drastically affect the “optimal” wing sizes and shapes taken on 
during the course of the design mission. This substantiates the need for a credible wing weight prediction early in 
the conceptual design phase to establish the necessary morphing capabilities. 

The conceptual idea behind the development of the morphing wing weight equation is relatively straightforward. 
However, implementation of the overall process unveils several problematic areas that must be approached with 
rigor to establish credible wing weight estimates. For example, the structural sizing process (the simultaneous 
analysis optimization approach) must be implemented to avoid the potential over-prediction caused by sequential 
optimization approaches. The method to perform this process is in place for the intermediate complexity model but 
has not yet been implemented. The validation and feasibility of this approach is of paramount importance and will 
be addressed as soon as an appropriate FEA environment is in place. Of further concern is the modeling of 
morphing wing skins in the FEM (skin material that undergoes strain as a direct result of rearrangement of the 
substructure). Currently, this issue is not formally addressed as it is unclear how to model the material behavior, 
load carrying capabilities, and failure modes. Although not addressed here, these materials are being developed, 
tested, and, in the MAS program, incorporated into morphing structures. For example, shape memory polymers 
were used in Lockheed’s folding wing design 3 ; the NextGen concept incorporates a “flexible skin for the wing that 
. . . smoothly accommodate(s) desired changes in the surface area 20 ”; and a recently-initiated DARPA program will 
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research and develop materials having characteristics similar to the cellular structure and behavior of plants (nastic 
materials) for use in adaptive / morphing structures. 21 Finally, the actuation system has been modeled here using 
simplified hydraulic actuators. Although this type of actuation scheme is conventional, a morphing wing would 
more likely use compact hybrid actuators 22 ' 23 as a means of weight and volume savings to increase overall system 
feasibility; these actuators will have increased power density and reduced volume compared to conventional 
actuation systems. However, weight estimates of these systems are currently unavailable or proprietary and can not 
be incorporated in an academic environment. Furthermore, this work assumes that the wing morphing is a quasi- 
static process with linearly varying loads. If the wing is desired to operate / morph within a high frequency range, 
the actuator sizing could be driven more by inertia and unsteady aerodynamic effects rather than loads due to 
friction and structural rearrangement. These concerns should be addressed in subsequent work. 

The template-based aircraft sizing routine enables sizing of an aircraft utilizing morphing wings with a specified 
mechanization strategy. This is important because resulting morphing wing capabilities will be physically 
realizable. An important consideration, however, is that the aircraft sizing routine requires on the order of one hour 
dedicated runtime (good initial guess for the design variables) to five hours (poor initial guess or over-constrained 
problem statement); these values are based on a dedicated desktop computer, Pentium IV processor (2.8 GHz), 
running Windows XP and Matlab 7.1 (R14). Techniques are currently being considered to expedite runtime 
performance by systematically decomposing the problem into a hierarchy of sizing-related subproblems; the overall 
goal and template-based modeling approach would be preserved, however. Finally, further studies are suggested to 
verify global optimality of the aircraft sizing results. While the obtained results seem reasonable, the optimality has 
not yet been established. 
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Appendices 


A. Examples of Template-Based Morphing FEM’s 


Rigid Body Variable Sweep Template 



Shear Sweep Template 
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Folding Wing Template 



B. Overall Process Implementation for the Beam Model Representation 


Environment Process 


Details 



DOE Type 

Matlab Command 

FC-CCD 

ccdesign(n, 'center', 1, 'fraction', 0, 'type', 'faced') 

FFD 

fullfact() 

D-opt 

rowexch(n, num_runs, 'purequadratic') 


Optimization algorithm : fmincon() 

FEA: In-house code 


Regression algorithm : IsqcurvefitQ 


Figure 9 Overall Process Implementation for the Beam Model Representation 
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C. Overall Process Implementation for the Intermediate Complexity Model 


iSIGHT 






DOE exp, i 

DOE 




i = i + 1 

i 


Morphing Wing 
FEM’s 


Simultaneous 
Analysis Structural 
Optimization 

1 


FEM decks with updated 
design variable values 




1 l 

Optimization 



ASTROS 

Algorithm 




Output File : 

• Objective Function 

• Constraint values 

• Strength 

• Buckling 

• Flutter Information 

• Gradient Information 


Optimal {a*} 



Figure 10 Overall Process Implementation for the Intermediate Complexity Model 
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